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As modified gravity theories, the 4-dimensional metric fCR) theories are cast into connection dynamical 
formalism with real .sj((2)-connections as configuration variables. This formalism enables us to extend the 
non-perturbative loop quantization scheme of general relativity to any metric fCR) theories. The quan- 
tum kinematical framework of fCR) gravity is rigorously constructed, where the quantum dynamics can 
be launched. Both Hamiltonian constraint operator and master constraint operator for /(K) theories are 
well defined. Our results show that the non-perturbative quantization procedure of loop quantum gravity are 
valid not only for general relativity but also for a rather general class of 4-dimensional metric theories of gravity. 

PACS numbers: 04.60.Pp, 04.50.Kd, 04.20.Fy 



O 

(N 

Oh! 

oo 



o : 



> 

(N 

On 
^, 

l> 

O 



X 



I. INTRODUCTION 



The theoretical search for a quantum theory of gravity has 
been rather active. Especially, as a background indepen- 
dent approach to quantize general relativity(GR), loop quan- 
tum gravity(LQG), has been widely investigated in recent 
twenty-five years. For reviews in this field, we refer toyj- 
0]. It is remarkable that, as a non-renormalizable theory, 
GR can be non-perturbatively quantized by the loop quanti- 
zation procedurep]. This background-independent quantiza- 
tion relies on the key observation that classical GR can be 
cast into the connection dynamical formalism with structure 
group of S U(2). Thus one is naturally led to ask whether GR 
is a unique relativistic theory of gravity with connection dy- 
namical character. Recently modified gravity theories have 
received increasinged attention in issues related to "dark en- 
ergy" and non-trivial tests on gravity beyond GR. A series of 
independent observations, including type la supernova, weak 
lens, cosmic microwave background anisotropy, baryon oscil- 
lation, etc, implied that our universe is currently undergoing a 
period of accelerated expansionfo]. This result conflicts with 
the prediction of GR and has carried the "dark energy" prob- 
lem. Although the acceleration could be explained by intro- 
ducing a cosmological constant A, the observed value of A is 
unnaturally much smaller than any estimation by tens of or- 
ders. Hence it is reasonable to consider the possibility that 
GR is not a valid theory of gravity on a cosmological scale. 
Since it was found that a small modification of the Einstein- 
Hilbert action by adding an inverse power term of curvature 
scalar K would lead to current acceleration of our universe, 
a large variety of models of fCR) modified gravity have been 
proposedfY]. Moreover, some models of fCR) gravity may 
account for the "dark matter" problem, which was revealed 
by the observed rotation curve of galaxy clusters. We refer to 
OTI ^ for a recent review on fCR) theories of gravity and It's 
application to cosmology. It is also worth noting that certain 
effective equation of loop quantum cosmology can be derived 
from some classical fCR) theory 1 9]. 
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Historically, Einstein's GR is the simplest relativistic theory 
of gravity with coiTect Newtonian limit. It is worth pursuing 
all alternatives, which provide a high chance to new physics. 
Recall that the precession of Mercury's orbit was at first at- 
tributed to some unobserved planet orbiting inside Mercury's 
orbit, but was actually explained only after the passage from 
Newtonian gravity to GR. Given the strong motivation to fCR) 
gravity, it is desirable to study such kind of theories at funda- 
mental quantum level. For metric fCR) theories, gravity is 
still geometry as in GR. The differences between them are 
just reflected in dynamical equations. Hence, a background- 
independent and non-perturbative quantization for fCR) grav- 
ity is preferable. The framework of extending LQG to fCR) 
theories appeared in 1 10]. The purpose of this paper is to pro- 
vide the detailed constructions. 

We will show that the connection dynamical formulation 
of fCR) gravity can be derived by canonical transformations 
from it's geometrical dynamics. The latter was realized by in- 
troducing a non-minimally coupled scalar field to replace the 
original fCR) action and doing Hamiltonian analysis. While 
the equivalence by canonical transformations at the classical 
level does not imply equivalence after quantization, our choice 
of the canonical formalism enables us to caiTy out the physical 
and mathematical ideas of LQG. The canonical variables of 
our Hamiltonian formalism of fCR) gravity consist of su{2)- 
connection A'^ and it's conjugate momentum E" , as well as 
the scalar field (p and it's momentum n. The Gaussian, diff'eo- 
morphism and Hamiltonian constraints are also obtained, and 
they comprise a first-class system. Loop quantization proce- 
dure is then naturally employed to quantize /(??) gravity. The 
rigorous Kinematical Hilbert space structure of loop quantum 
GR is extended to loop quantum /(!??) gravity by adding a 
polymer-like quantum scalar field. The spatial geometric op- 
erators of LQG, such as the area and volume operators are still 
valid here. Hence the important physical result that both the 
area and the volume are discrete at quantum kinematical level 
is also true for /(■??) gravity. As in LQG, the Gaussian and dif- 
feomorphism constraints can be solved at quantum level, and 
both the Hamiltonian constraint and the master constraint can 
be promoted to well-defined operators. 

This paper is organized as follows. In section|lll we derive 
the connection dynamical formalism for /(■??) theories. In sec- 
tion |III1 the kinematical Hilbert space for/('R) gravity is con- 



structed, where the Gaussian and diffeomorphism constraints 
are implemented. The Hamihonian constraint is promoted to a 
well-defined operator in the kinematical Hilbert space in sec- 
tion |IV] We also define a self-adjoint master constraint oper- 
ator in the diffeomorphism invariant Hilbert space in section 
rvl Finally, some concluding remarks are given in section [VTl 
We use Greek alphabet for spacetime indices. Latin alphabet 
a,b,c,...,for spatial indices, and i,j,k,..., for internal indices. 



II. CONNECTION DYNAMICAL FORMALISM FOR f(R) 
THEORY 

A simple extension of GR is to consider the Lagrangian 
of gravity as a function of scalar curvature fi. This kind of 
modified gravity theories have become topical in cosmology 
and astro-physics. The original action of fCR) theories read: 



S(g)^^-jd\^P§f('R) 



(1) 



where / is a general function ofR, and we set SnG = 1 . By in- 
troducing an independent variable s and a Lagrange multiplier 
(/), an equivalent action is proposed as llllL[l2ll : 

S(g,<P,s) = ^Jd'x^/^(f(s)-<p(s-'R))■ (2) 
The variation of dU with respect to s yields 



ds 



^ /'(*)• 



(3) 



Assuming f"(s) ^ so that s could be resolved from the 
above equation, action ^ is reduced to 



S(g,<P) = ^jd\^f^(<p'R-^m = Jd'x£(. 



X) (4) 



where ^((p) = <ps - fis). The variations of (HJl give the follow- 
ing equations of motion 



n = r(0) 



(6) 



where ^'((f>) = -§P-, and V^ is the connection compatible with 
gfjy. It is easy to see that Eqs. (|5]l and (|6]l are equivalent to 
the equations of motion derived from action ([T]l- The virtue of 
action (|4]i is that it admit a treatable Hamiltonian analvsis Lllll . 
The Hamiltonian formalism can be derived by doing 3-1-1 de- 
composition and Legendre transformation: 



^ab 



dZ 

d'hab 

^ WK"' - Kh"') - ^(0 - A^^5c0)], (7) 



Uah 



d(j) 



(8) 



where hah and Kat, are respectively the induced 3-metric and 
the extrinsic curvature of the spatial hypersurface E and K = 
K^. The combination of the trace of Eq.(|7]i and Eq.® yields 



2N 
■N'dccf>^^(<f>n-p). 

3V/! 



(9) 



Note that the action (01 is also meaningful for a constant (p. 
In this special case, one could resolve n from p by Eq.® as 
n = p/if). This reduces one degree of freedom of the theory. 
Then the fCR) theory will be reduced back to GR. In general 
case, the Hamiltonian of fCR) gravity can be derived as a liner 
combination of constraints as 



Hfotal - I <: 



d'x(N"Va+NH). 



(10) 



where A^ and A^" are the lapse function and shift vector re- 
spectively, and the smeared diffeomorphism and Hamiltonian 
constraints read 



V(N) = I d\N''Va 






d'xN"(-2D''ipab)+nda<l>), 



(11) 



H(N) = d'xNH 

r -, 2 pabp"'' -W 1 , 1 

J-z ^Th (J) 6 3 

+ ^V^(^(0)-0/? + 2D„D»], (12) 

where Dg is the connection compatible with the 3-metric hat,- 
Again, in the special case of <;* = costant, it is easy to see that 
the smeared diffeomorphism and Hamiltonian constraints can 
go back to GR up to a constant rescale. By the symplectic 
structure 



r(crrf)c3/ 



{K,h{x),p'%y)} = ^^^,'6\x,y), 
{<P(x)My)] = ^'(^,3'), 



(13) 



lengthy but straightforward calculations show that the con- 
straints (fTTT i and (fT2l i comprise a first class system similar to 
GRas: 



{y(A^), y(A^')) = y([A^,A^']), 
{V(n\H{M)) = mx^M), 

{H{N),H{M)} = V{ND"M-MD"N). 



(14) 



Since the above Hamiltonian analysis is started with the ac- 
tion dU where a non-minimally coupled scalar field is intro- 
duced, we need to check whether the Hamiltonian formalism 
is equivalent to the Lagrangian formalism. It is not difficult to 
see from the Hamiltonian (fTOl) that the evolution equation of 
the scalar field reads: 



{cl>,H, 



total I 



2N 
3^fh 



(<Pn-p) + N"da 



(15) 



which is nothing but the Eq. of (|9]). The evolution equation of 
the 3-metric reads 



A^ Mpab - \phab) 2 
hah = —p{ -r -nh„b) + DaNh + DbNa 



INKah + DaNb + DbNa, 



(16) 



which is nothing but the definition of Kab ■ The evolution equa- 
tion of the momentum of reads 

2N ^ PabP"'' -\p^ 1 _2, N^^, 



-^n^)-—m 



n = 5„(A^V) + — ( 
yfh 

NVh rr 

+ — —R - ^iDaD"N 

2 

+ E2^R-da{ylhh''''dbN). 



(17) 



Using the definition of tt = - ^K and rP - jj,n' - -^, we 
can get 



r(0) = KabK"" - K^ + R + -d.i^f^n'K) 

N^/h 



N^fh 



(18) 



This is nothing but Eq. (|6]l. On the other hand, the 00- 
component of Eq.© reads 

cl>G,Xrf = '^{R - KabK"" + K^). (19) 

Using the identity g^yrfrf — -\, Eq.(fT9]l becomes 

'^{R-K„bK"' + K') 

= ^^(0) + (g^'' + «'^nOV^V,^ 

= ^^(0) + D„D>-/:(i(0-A?^5.<^)) (20) 

where the facts h^'^Uy - and rfda-(p - j^{^ - N'^dcCp) have 
been used in the above derivation. Note that the Hamiltonian 
constraint in Eq.(fT2t can be expressed as 

H = 2El(K,bK'"-K^-R) + ^(^(<P) + 2Daiycl>) 



VhKi^(i>-N'd,<p)). 



(21) 



Hence the 00-component of (|5]l is equivalent to Hamiltonian 
constraint. Now we come to the Oa-component of (|5]). Since 

4>G,yn''hl = cp{DaK-D,K) and g^^n'^hl^Q, (22) 

we have 

m^Kl - DhK) = nlilV^V,4 

= «X'^cr((/i^ - «V)5^<^). (23) 



The first term in the right hand side of above equation reads 



nli'^WAh'id^'l') = nli:WA/y+n^ny)d^4 



-K'A 



a"b(l 



(24) 



and the second term reads 



n^hZVJn^'ny^,^) = h'^^^Jn'^d,^) 



= Da{-(4>-N'dc4>)). (25) 



Hence their combination gives 



DMKD - Db(cf>K) + KDbcf) - Db(-(cf> - N'd, 

2 y[h K . 

^Da( — [4>(Kl - Khl) - -^(c/, - N'dcc/>)]) 

yh ^ ^ 



yfh 



dhcp. 



(26) 



This is nothing but the difi'eomorphism constraint in Eq.lfTTTi. 
Now we turn to the ab-components of Q. We will show that 
they are equivalent to the equation of motion of pab which 
reads 

habN PcdP"" - \p^ 1. 1 
Pab = — ^( 7 + t4>^ - oP^) 

yfh fp 6 3 

2N PacPl-\pPab 1 
+ -^( 7 -;Pab7^) 

+ - ^fhhab(l>R - J ^(pRab - J ^Kb^icP) 

- 2" ^habD.D'cp - D(„N ^Dh)(P 

yfh 
+ —(D(,Db)(N<p) - habD,D'(N<p)) 

+ 2p,(aD'Nb) + D,(PabN% (27) 

Since the initial value formalism of /("??) theories has been ob- 
tained in fTJl, we will use Eq.jZTTi to derive the time derivative 
of the extrinsic curvature: 



Kab - 



^(Pab - \phab) nh 



ab 



(28) 



ab 



4>yfh 3V/i 

A straightforward calculation yields 

kab = 2NKacKl-NKK„b + J:^K„b-NR. 

+ DaDbN + -D„Dhcf> 
<f> 

Nhab , f(0) n'^dir<p 

+ ^(^'(<f>)+^-^)-—f^NKab. (29) 

It is easy to see that Eq.(|29]l is equivalent to Eq.(217) in Q. 
Note that there is a sign difference between the definition of 
our extrinsic curvature and that in |7], and our potential ^(0) 
is twice of that in [7]. To summarize, we have shown that the 
Hamiltonian formalism of fCR) gravity is equivalent to it's 
Lagrangian formalism. 



Recall that the non-perturbative loop quantization of GR 
was based on it's connection dynamic formalism. It is very 
interesting to study whether the previous geometric dynamics 
of fiH) modified gravity also has a connection dynamic cor- 
respondence. To this aim, we first extend the phase space of 
geometrical dynamics to the triad formalism, and then intro- 
duce a canonical transformation on the extended phase space 
of/(!R) theories. Let 



Uab 

IN 



= <t>K"'' + 



uab 



3^fh 



t>n - p). 



(30) 



we obtain the standard Gaussian constraint 



Qi = ^El = daEl + eijkAiE 



ak 



(37) 



which justifies A'^ as an iM(2)-connection. Note that, had we 
let y - +i, the (anti-)self-dual complex connection formalism 
would be obtained. The original diffeomorphism constraint 
can be expressed in terms of new variables up to Gaussian 
constraint as 



Va = -2D''ipab) + ndc4 
- \p'abEl+^d,'P, 



(38) 



and £■? = yhe", where e? is the triad s.t. hahs'le''. - Sij. Then 
we get 



,^ab 



= ^ikfE"' - ^-K[E^,E)E% 



(31) 



where K" = K^^'e'i^. Now we extend the phase space of geom- 
etry to the space consisting of pairs (E", K'^). It is then easy to 
see that the symplectic structure (fTsT l can be derived from the 
following Poisson brackets: 



{E"j(x),E',(y)} = {Kiix),Kl(y)} = 0, 
{Ki(x),E',(y)}^d'„5id(x,y). 



(32) 



Thus there is a symplectic reduction from the extended phase 
space to the original one, and the transformation from con- 
jugate pairs (hab,p"') to (£", A'^) is "canonical" in this sense. 
Note that since K"'' - &", we have an additional constraint: 



Gjk = ka\jEl^ - 0. 



(33) 



So we can further make a canonical transformation by defin- 
ing: 



Ai = rj, + r^. 



(34) 



where FJ, is the spin connection determined by E", and y is a 
nonzero real number. It is clear that our new variable A'^ co- 
incides with the Ashtekar-Barbero connection lll3Lll4ll when 
0=1. The Poisson brackets among the new variables read: 



{Ai(x),E'liy)] = y5'X5{x,y), 
{A[{x),Ai{y)] = 0. 



(35) 



Now, the phase space of /(■??) gravity consists of conjugate 
pairs (AJ,,/?^) and (0, tt). Combining Eq.(l33Tl with the com- 
patibility condition: 



daE^; + eijkViE"' = 0, 



(36) 



where F^^ = 2d[aA' + e'l^fA^Ai^ is the curvature of AJ,. The 
original Hamiltonian constraint can be written up to Gaussian 
constraint as 



H = i;iKb-(r + i-2')^M'KK] 



S'ln i>n -I J" k 



1 2 (KEff 4 {kiE^)n 2 n^ 
+ -( ■ 1 ■ 1 

2 ^<P yfh ^ yfh ^ yfh 

+ V/^(0)) + VhDaD''(p. 



(39) 



It is easy to check that the smeared Gaussian constraint, 
^(A) := J d^xA'(x)Gi(x), generates SU(2) gauge transfor- 
mations on the phase space, while the smeared constraint 



'V(N) 



-L 



d\N''(V„-A[gd 



(40) 



generates spatial diffeomorphism transformations on the 
phase space. Together with the smeared Hamiltonian con- 
straint H{N) - J cPxNH, we can show that the constraints 
algebra has the following form: 

{^(A),^(A')) = ^([A,A']), (41) 

{g{K),^(N)] = -gU^K), (42) 

{g{K),H{N)} = 0, (43) 

{^{N),^(n')] = ^{[N,'n']), (44) 

{^(N),H(M)] = HiL^M), (45) 

{H(N),H{M)) = 'V{ND"M-MD"N) 

+ g{(Nd„M-Md,N)h"'Ah) 
[E''D,N,E''DbM]' 

- -, gi 



2[E"D„(<f>N),E'Dh((f>M)Y 



7 



gi. (46) 



One may understand Eqs. ( 14 11143] ) by the geometrical interpre- 

— > 
tations of ^(A) and ^(A^). The detail calculation on the Pois- 
son bracket ( l46b between the two smeared Hamiltonian con- 
straints will be presented in the Appendix. Hence the con- 
straints are of first class. Moreover, the constraint algebra of 
GR can be recovered for the special case when (p - I. The 



total Hamiltonian is a linear combination of the above con- 
straints as 



n. 



L 



H(N)+N"V„+g(A). 



(47) 



To summarize, fCR) theories of gravity have been cast 
into the iM(2)-connection dynamical formahsm. Though a 
scalar field is non-minimally coupled, the resulted Hamilto- 
nian structure is similar to GR. Note that what we obtain 
is real sM(2)-connection dynamics of Lorentian fCR) grav- 
ity rather than the connection dynamics of some conformal 
theoriesHQ. 



III. QUANTUM KINEMATIC OF fCR) THEORY 

Recall that LQG is based on the connection dynamics of 
GR. We have shown in last section that fCR) theories can also 
be reformulated as connection dynamical theories by intro- 
ducing a non-minimally coupled scalar field. Hence the non- 
perturbative loop quantization procedure can be straightfor- 
wardly generalized to fCR) theories. Since the configuration 
space consists of geometry sector and scalar sector, we expect 
the kinematical Hilbert space of the system to be an direct 
product of the Hilbert space of geometry and that of scalar 
field. To construct quantum kinematics for geometry as in 
LQG, we have to extend the space £/ of smooth connections 
to space £/ of distributional connections. A simple element 
A e £/ may be thought as a holonomy. 



hJA) = !Pexp 



I'" 



(48) 



of a connection along an edge e c 2. Through projective 
techniques, ^ is equipped with a natural measure fiQ, called 
the Ashtekar-Lewandowski measureylQ]. On the other hand, 
one may smear the densitied triad E" on 2-surfaces to obtain 
fluxes as 



E{S,f) 



'■- I eabcE"f' 



(49) 



where /' is a ,s'M(2)-valued function on S . From the algebraic 
viewpoint, the cylindrical functions of holonomies and the 
fluxes consist of an C*-algebra. Then by Gel'fand-Naimark- 
Segal(GNS) structure|2], one can obtain the cyclic represen- 
tation for the quantum holonomy-flux *-algebra, which co- 
incides with the one by projective techniques. In a cer- 
tain sense, this is the unique diffeomorphism and internal 
gauge invariant representation for the quantum holonomy-flux 
algebral 17]. The kinematical Hilbert space of geometry then 
reads l-l^^^ - l}(s^ , dii()). A typical vector ^a{A) € 'TYy^ is 
a cylindrical function over some finite graph a c S. The so- 
called spin-network basis 



Ta(A) = Y[ V2}rn'<,„/A(e)), ije^Q) 



eeE(a) 



provides an orthonormal basis for "T^yjilj, 0], where 
tt/,, „ (A(e)) denotes the matrix elements in the spin-y repre- 



sentation of S U{2). Note that the spatial geometric opera- 
tor of LO G, such as the areajlSJ] , the volume|19] and the 
lengthl20l ulll operators, are still valid in "H^'^, though their 
properties in the physical Hilbert space still need to be clari- 
fied 1 22, 23]. 

Since the scalar field also reflects fCR) gravity, it is natural 
to ernploy the polymer-like representation for it's quantiza- 
tion [24]. In this representation, one extends the space ^ of 
smooth scalar fields to the quantum configuration space ^. A 
simple element U ^'% may be thought as a point holonomy. 



Ux - exp(//10(x)). 



(51) 



at point X € S, where /I is a real number By GNS structure 113], 
there is a natural diffeomorphism invariant measure dyi on 
'% I24i] . Thus the kinematical Hilbert space of scalar field 
reads 'H^^ - L?{'^,dfi). The following scalar-network func- 
tion of (f): 



Txiip) = Tx,M) = n ^-i^^^^^))' 



(52) 



x,eK 



where Z = {xi, . . . , x„) is an arbitrary given set of finite num- 
ber of points in S, constitute a complete set of orthonormal 
basis in ^^^^- Since the point holonomy of a scalar is defined 
on an 0-dimensional point, the momentum is smeared on an 
3 -dimensional region /? in Z as; 



n{R) 



-L 



d xn{x). 



(53) 



Thus the total kinematical Hilbert space for /(■??) gravity reads 
"Hkin := "H^^ ® TY^Hj with an orthonormal basis Ta,x{A,(p) = 
Ta{A) ® Tx(4>)- Note that a basic feature of loop quantization 
is that only holonomies will become configuration operators, 
rather than the classical configuration variables themselves. 
Let ^(A,^) denote a quantum state in '74in. The actions of 
basic operators read 



h,(AW(A,(f>) 
E(S,fmA,cf>) 

umv)mA,4>) 

7i(R)'i'(A, (f,) 



hMmA,4>), 

in{E(S,f),nA,cf>)}, 
exp(iA(f>{v))^'(A, (p), 
in{n{R),'V(A,4))]. 



(54) 



(50) 



As in LQG, it is straight-forward to promote the Gaussian con- 
straint ^(A) to a well-defined operator in 'Kidn. It's kernel is 
the internal gauge invariant Hilbert space 'Kg with gauge in- 
variant spin-scalar-network basis Ts,c - Ts{A) ® Tx((p), where 

Ts=(a,jdA) = ®vev(a)iv ' m^'{A{e)), (;, ?t 0). (55) 

Here an intertwiner / is assigned to each vertex of graph a. All 
the internal gauge invariant geometric operators, such as the 
area, volume and length, can also be well defined in 'He- Since 
the diffeomorphisms of Z act covariantly on the cylindrical 
functions in "He, the so-called group averaging technique can 
be employed to solve the diffeomorphism constraintt3n4J. To 



this aim, we first define a projection map acting on cylindrical 
functions if/p = tlf^^xiA, (p) e "Kkin as 



1 Y^ . 



(56) 



where t/^ denotes the unitary operator corresponding to a fi- 
nite diffeomorphism i/j : 2 — > S. Here GSp - DiffplTDiffp 
is the group of graph symmetries, where Diffp is the group of 
all diffeomorphisms preserving the colored y6, TDiffp is the 
group of diffeomorphisms which trivially acts on fi, and tip 
is the number of the elements in GS p. Secondly, we average 
with respect to all remaining diffeomorphisms which change 
the graph (3. For each cylindrical function ilip, there is an el- 
ement rjiij/p) associated to it in the algebraic dual space Cyl* 
which acts on any cylindrical function (pp' as 



ipeDiff(J.)IDiffy 



(57) 



Jl 30 ^|h 

Ji 3 V^ 

r . Nn'^d) 

Js 3 V^ 



(60) 



Here we have written the smeared version of Eq.([39]l as 
H{N) - Tj]=\Hi. Note that the first two terms in H(N) can 
be written as 



Hi - ^- j^d^xNcpF^, 



'3„Arj,ri ^Jl''^k^l 



V/^ 



H'^"''((/)N), 



(61) 



where DiffiY) is the diffeomorphism group of 2. It is easy 
to verify that ri{\jjp) is invariant under the group action of 
DiffCL), since 



vW[U^(f 



rji^PfiMfi']- 



(58) 



Thus we have defined a rigging map rj : Cyl — > Cyl'^.r,, which 
maps every cylindrical function to a diffeomorphism invariant 
one. Moreover, a Hermitian inner product can be defined on 
^y^Diff ^^^ ^^^ natural action of the algebraic functional: 



(nW\v(^/S'))Diff := rii^/sM/}']- 



(59) 



The diffeomorphism invariant Hilbert space 'Hoiff is defined 
by the completion of Cyl'^..^ with respect to the above inner 
product. Thus we can also obtain the desired diffeomorphism 
and gauge invariant Hilbert space, ^oiff, for /(^) gravity. 



Hi 



--IL 



1 eiiE"E'' 

1 ._ 1 



1+y^ 



-T{N{j'^ + -)), 



(62) 



where T{N) denotes the Lorentzian term in the Hamiltonian 
constraint of GR. Hence, except that the smearing functions 
are multiplied by some function of (p, these terms keep the 
same forms as those in GR. 

By introducing certain small constant Aq, an operator cor- 
responding to the scalar (p{x) at x e S can be defined as 



ZIAq 



(63) 



The ambiguity of Aq is the price that we have to pay in order to 
represent field in the polymer-like representation. To further 
define an operator corresponding to (p^^{x), we can use the 
classical identity 



IV. QUANTUM HAMILTONIAN OF fCR) THEORY 

While the kinematical frameworks of LQG and polymer- 
like scalar field have been straight-forwardly extended to fCR) 
theories, the nontrivial task is to implement the Hamiltonian 
constraint ( [39] l at quantum level. In this section, we can show 
by detail and technical analysis that, as in LQG, the Hamilto- 
nian constraint can be promoted to a well-defined operator in 
the kinematical Hilbert space "Kkin- The resulted Hamiltonian 
constraint operator is internal gauge invariant and diffeomor- 
phism covariant. Hence it is at least well defined in the gauge 
invariant Hilbert space "He- 

Comparing Eq.([39]l with the Hamiltonian constraint of GR 
in connection formalism, the new ingredient of fCR) gravity 
that we have to deal with are 4>(x), (p^\x), ^((p) and the follow- 
ing four terms 



-'(x) = sgnWijSgnWm'ixlniR)})- 



(64) 



for any rational number I 6 (0, 1), where sgn[(/i] denotes the 
sign function of 0, \<p\ is the absolute value of and x e R. For 
example, one may choose I - j for positive 0(x) and replace 
the Poisson bracket by commutator to define 



-'(x)^(^[^4,(xmR)]f 



(65) 



Thus all the functions ^((p) which can be expanded as pow- 
ers of 4>(x) have been quantized. For other non-trivial types 
of ^(0), we may replace the argument (;* by in Eq.(l63]l. pro- 
vided that no divergence would arise after the replacement. 
In the case where divergence does appear, there remain the 
possibilities to employ tricks similar to Eq.(l64ll to deal with 



it. Hence it is reasonable to believe that most physically in- 
teresting functions ^(0) can be quantized. Then it is straight- 
forward to quantize H(, = j Ln yh ■ ^(^) as an operator acting 
on an basis vector Tax as 

H6 ■ T,,x = ^ 2 ^(^')^(^(^))^>' ■ Ta,x- (66) 

veV{a) 

Note that the action of the volume operator V on a spin- 
network basis vector T^iA) over a graph a can be factorized 

as 



V-T, 



'a = X ^•' ■ ^«- 



(67) 



veV(a) 



Moreover, by the regularization techniques developed for the 
Hamiltonian constraint operators of LQG and polymer-like 
scalar field, all the terms Hi,H4,Hs and H-i can be regular- 
ized as operators acting on cylindrical functions in 'Kkin in 
state-dependent ways. In the regularization procedure, we 
will use the following classical identities 
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(68) 



where K - L cPxK'^E" can be write as Poisson bracket: 



^ = y-I{H^"'''(l),y). 



(69) 



Here the Euclidean scalar constraint H (\) by definition 

was: 



//^'"^'(l) 



- \\/<^ 



^jmEIEI 



yfh 



(70) 



Both H^"'^' and the volume V under consideration have been 
quantized in LQG. Also, one has E" - 5eyte"'"'e^e^, where 
e"*'^ is the levi-civita tensor density, and the co-triad satisfies 
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(71) 



To deal with the four new terms ( |60b , we first regularize them 
separately by point-splitting and obtain 



Hi = Hm 



sX-'-'i 



4 A^ 



d x—-x^(x - y) 



K',(x)Ef(x) Kl{y)E]{y) 






{H^"''il), .JiV^)}, 



(72) 



H4 = lim 



i!i r 



34/ 



d y7:(y)xAw - y) 



I 



d'xNx.(x-y){A'Jx),K} 



Mbc 



X e'"'^Tr(T,{A,(x), (Vvir''][Ac{x), (Vc/j)^'")) 



^3/4> 



I' 



Jefr, 






r , C -x N(h 

Hs ^ hm \ d y \ d x—-xJx - y)7T(x)n(y) 

r 3 det{e'^(u)) 
X I (i'^M- 



(^^„)'^' 



-;^e(m - x) 



3 flfef(e;,(w)) 
fl w ,.., .,,. ;\f^(w - y) 



= lim 



Je Ji 



d xN(pTi{x)n{y) 



. -0 3^^ 

^^^(^ - y)^^^ (m - -^)a'£("^ - y) 



I li^Mi 



_a/7c 



X I d"Me°''"Tr({A„(M), VVf/j) 
{Ai(M), Vv^){A,(m), Vv^)) 

X rfl!Ve''^^Tr({Arf(w), V^^) 
{Ae(w), V^){A/(w), V^)), 



?£>/,0 



(74) 



I { d^y { d\{D^,(x-y)) 

' J-E. J-E. 



- lim 



X N{x)E'^{x) 
X (Di,(l>(y)) 



Efiy) 



e3(/i(3;))l/2 



-lim 



r^ J2 J2 



fl'^x(D„;i^,(x-y))A?(x)£f(x) 



X e''HDbcP(y))TT(TdAAy), (Vf/,0'^'){A/(y), (Vt/.O'^'l), 

(75) 

where ;t'£(x - y) is the characteristic function of a box {/^ con- 
taining X with scale e and satisfies the relation limg^oXeix - 
y)le^ - 6{x - y), and Vu^ denote the volume of U\. It is easy 
to see that the regulator in i/3 can be removed by acting on a 
given basis vector Ta,x e "Wkin as 



^3 ■ Ta,x 



AN{v) . 1 



'ijHiny 



r'(v) 



2 

veV(a) 



X [//£-'(!), VK][//£'"'(1), VK] ■ r.,z, 



(76) 



For the other three terms, in order to reexpress connection 
by holonomy and make the regularization diffeomorphism co- 
variant, we triangulate E in adaptation to some graph a under- 
ling a cylindrical function in "Hkin. At every vertex v e Via), 
for each triple (e/, ej, ck) of edges of a we have a tetrahedron 
Af;„ „ „ basedatv, which is spanned by segments ,?/, ,sy, i/f of 
the triple. Each segment si is given by the part with the curve 
parameter t' e [0, e] of the corresponding edge ei(t'). For 
each A*; „ _ „ one can construct seven additional tetrahedron 
(73) by backward analytic extension of the segments. The regions 



of Z without a vertex of a can be triangulated arbitrarily. Note 
that for one segment sj, we have 



f{A(u), ^JV(u,e)} ^ ei^(0){A„(v), ^JV(u,e)} (77) 

U Sf 



respectively 
i/^ = -lim 

6— 

X Z 

l'GQ'(l') 



i20 



X 



3*y* 

A^(v(A)) 
£(v) 



d yniy) 



up to 0(e). Hence for each AJ^ , 



, we have 



i'"e»(v) ^ ' v(A)=v" 



X /z 



r e^^^rr ({A„(m), VW } [Ab(u), ^|vf, } {a,(m), VW }) 






//f = lim- 



(78) 



-0 35y 



5 ,,6 



X 



I £/^xA?(x)0(x)7r(x) I d^^ 



y7r(y) 



^2 JS 

X ;r£(v(A"') - y);(^^(y(A") - x)xeix - y) 



where eisjSjSK) '■- sgn(det(i/iyijs-)(y)) takes the values 
+ 1,-1,0 if the tangents of the three segments si, sj, sk at v 
(in that sequence) form a matrix of positive, negative or van- 
ishing determinant. Then the integration over E can be split as 
follows ill: 



ZeiSrS rSj()e' 

v"ea{v) ^ ' v(A)=v" 



Z 



UK 



X Tr(/t 



■'/(A")f''.J(A")' -\ 



Vu%. ]h 



^XA"){''j,(A")' ^| 



Vi, 



X h. 



Ah 



■5x(A )l".5jf(A")^ 



v^» 



J- Zj e(v) ^ 

- { -y — y 

\,y, Zj e{v) Zj 



X 



z 



1 



E(v"') Z e(^L^Miiv)e^* 

i'"'eff(i') ^ ' i<A)=r"' 



i(«,)ni(o)ni)(eK)=i' L 



veyCff) 



b{e,)^\b(eJ)(^b{eK)=v 






V^)/^ 



X Tr(/i,^(A"'){/!;^(^,.), 



Ah 



-1 



•5m(A )l' j„(A"') 



'V^ 



(79)//,- = -lim 



l''Eff(lO '^ ' v(A')=l'' 



fd'x(DaXe(x-v'))N(x)Ef(x) 



SjSk 



•^^JK 



Here we have first decomposed E into a region {7^^ not con- 
taining the vertices of a and the regions Uf-, ,, around the ver- 
tices. Then choose a triple (e/, ej, e^) of edges outgoing from 
V and decompose t/^ ,, into the region U^ ,, ^^ ^,^ ^^ covered by 
the tetrahedron AJ; „ „ „ „ spanned by ei,ej, bk and its 7 mir- 
ror images and the rest U^ ,, ^^ ^.^ ^,^ not containing v. Note that 
the integral over U^ ^ ^^ ^^ ^^ classically converges to 8 times 
the integral over the original single tetrahedron A^ ,, ^^ ^^ ^^ as 
we shrink the tetrahedron to zero. We average over all such 
tripes (e/, ej, ck) and divide by the number of possible choices 
of triples for a vertex v with n(v) edges, E{v) - ("3'). Then 
by the above triangulation T(e), the regulated 3 terms become 



X t/^„'(0(5,,(A')))[t/i„(0(f.,(A'))) - t/^(0(^.>/(A')))] 

X Tr(T,^,,(A,){/z;;(^,,, {Vui,f'^}hs,(A'){h-l^^,y {Vui,)"^}), 

(80) 

where v(A) and >s'/(A) denotes a vertex and a segment of a 
tetrahedron respectively, and f 5,(A) ( ■s.s,(a)) denotes target point 
(starting point) of a segment sjils). Note that the action of the 
operator n(R) on a scalar-network basis vector Tx{(t>) over a 
graph X can be factorized as 



n{R)-Tx^ 2 7i,rTx 



(81) 



.t,exn« 



Now every ingredient of H': has clearly quantum analogy, we 
can define the corresponding operators acting on a basis vector 
r„ V over some pranh nli Y as 



Ta^i over some graph a U X as 



Ha ■ Tax - - lim — r-^ r 

X ^ n{v')xe{v" - v')xe{v' - v) 

v'eX 

^ Z §S Z Tr(T,./..,(..,[/.;;,,,,,i]) 

veQ-(v) i'(A)=v 



X e 






UK 



v"eQ'(v) v(A)=v" 

X Tr(/!„(A„)[/i;/(^„),(yy.„)'/2] 



(82) 



2l8 

Hf . r„x = lim ^—^ — T y y 

X ^(v)N(v)7t(v)7t(v')x,(v"' - v')x,(v" - v)x,(v' - v) 



X 



E(v") 

v"ea{v) v(A)=v" 



X Tr(/i,,(A„)[/i;,U,„)>(%,)'^'] 

X ^.,,(A,.)[/^;,W,)'(^^:J'^'] 

X /j..(A.,)[^;V")'(^^,!.)'^']) 

X Z ^ Z K^.^M..)e^-- 

v"'eQ'(v) v(A)=v"' 

X Tr(/i,,(A,„)[^;/,^^^,^^),(K.O'^'] 

X /i.„(A,„)[/i;X,)' (\'„, )'^']) ■ Ta,x, (83) 



X IxAeitk) - v') -Ar.(e(ft-i) - v'Wieih^i)) 

X Z ^ Z -(^/^^^-)^'" 

v'eo-Ci') ^ ^ v(A)=v' 
X j7j^'(0(s.,(A„o))[^^o('^(f.v,(A,.,))) - UaM^MA,)))] 

X Tr(T,/i,,(A,)[/i-;^^,),(yy.,)'/2] 
-1 

.sk(A,/)' 



X /i.v^(A..)fc\A.),(V'c/',)'^'])-r„,z, 



(84) 



where = fo < fi < .. < f„ = 1 is an arbitrary partition of 
the interval [0, 1], Xi(t) := [hAO, t)TMt, l)]ABd/d[hAO, 1)]ab 
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(we denote X'^ := X^(0) in the following), and /zs,(a,.) denotes 
the holonomy along the segment sj starting from the vertex 
V of tetrahedron A. On the other hand, for H!j, we perform 
the limit n ^ oo, and e — > in reversed order Keeping n 
fixed, for small enough e, only the term with fe = 1 in the sum 
survives provided that i/(0) = v'. So for small enough e, the 
above operator reduces to 

m-Ta,x = lim 



. .0 3j^iAQ(ifi)^ 
X Y, ^^(0);r.(e(0) - v')N(em 



eeE(a) 



UK 



X Z ^ Z ^^^'^j^kW' 

v'eff(v) ^ ^ y(A)=v' 
X ^,l(0(s,,(A„,)))[£>i„(^(f.,(A„,))) - ^^o(<^(«.v,(A„,)))] 



X Tr(T,^,,(A,„)[/i-;(^^,),(yy.,)'/^] 



1/2-, 



X /l.v^(A..)[^;;(A^,,),(Vc/«,)"1)-r,,x. 



(85) 



Since the actions of H^ and H^ are independent of e, we can 
take the limits and obtain 



Ha ■ Taj 



E 



2^'W(v) 



3''/('S)''£^(v) 



Tt{v) 



veV{iy)nX ' 

X ^ Tr(T,/i,,(A)[/i;;(A),i]) 

v(A)=v(A')=i' 

X 6(,Vi.M^^)e^*''' 

X Tr(T,/l.„(A)[/i;^(A)- (K)'/']^.„(A)[^;J(A), (K)'^']) 

X e(.?/iyix)e^-^^ 

X /i.,,(A')[/i;i(A')'(^")'^'])-^«.^' (86) 



v^ 2"*A^(v) ^ . 
Hs-T„x — / ^ ^ 7 — ^ — f[(v)d)(v)n(v) 



X 2_] £(SjSjSk) 

v(A)=v(A')=i' 
X 



jyA: 



X /j..(A)[/i;iA),(K)'/']) 

X e(sLSMSN)e 



XMN 



X Tr(^.,,(A')[/i.,'(A') 



,,),(y,)^^']/^.„(Ao[/^;:(A')'(^")'^'] 

X hs„m[h-'^^,y(V.)"^])-Ta,x- (87) 

However, it is easy to see that the action of H^ on T^x is graph 
changing. It adds a finite number of vertices at t(sj(v)) = s for 
edges e/(f) starting from each high-valent vertex of a. As a 
result, the family of operators Hj(N) fails to be weakly con- 
vergent when e — > 0. However, due to the diffeomorphism 
covariant properties of the triangulation, the limit operator 
can be well defined via the so-called uniform Rovelli-Smolin 
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topology induced by diffeomorphism-invariant states (t>Diff 
as: 



ODiff(H^ ■ r„,x) = lim((Do,//|i/7^|r„,z). 



(88) 



It is obviously that the limit is independent of e. Hence both 
the regulators e and s can be removed. We then have 



Hi ■ Taj 



veV(a) 



3y^iAo(in)^E(v) 






Y^KY^ eisiSjSK)^' 



e(0)=i' i'(A)=y 
V-1/ 



X t/^;(</'(s.„(A)))[t/io(0(f.<,(A)))- i/^„(0(s..,(A)))] 
X Tr(T,-^,,(A))[/i;;(A),(K)'^'] 
X ^.^(A)[/2;^'(A),(K)'^'])-r,,x. 



(89) 



Collecting all terms, the whole Hamiltonian constraint can be 
quantized as a well-defined operator fKN) in 'Kkin- The action 
of li{N) on Tax can be factorized as 



^(A^) ■ Ta,x = Yj ^^^^'- ■ ^'^- 



(90) 



v£V(a) 



This operator is internal gauge invariant and hence also well 
defined in "Hg- However, although H(N) can dually act on the 
diffeomorphism invariant states, there is no guarantee for the 
resulted states to be still diffeomorphism invariant. 



V. MASTER CONSTRAINT OPERATOR 

Although the Hamiltonian constraint operator constructed 
in last section is well defined in 'Hq, it is diflicult to define 
it directly on "Hoiff- Moreover, the constraint algebra ( 1411 1- 
( |46] | do not form a Lie algebra. This might lead to quantum 
anomaly after quantization. In order to avoid possible quan- 
tum anomaly and find the physical Hilbert space, master con- 
straint programme was first introduced by Thiemann in IZSIl . 
We now apply this programme to quantum fCR) gravity. 

By definition the master constraint of fCR) theories classi- 
cally reads 



M 



■-'d 



,3 m^r 

a X- 



^/h 



(91) 



where the Hamiltonian constraint H(x) was given by Eq. 
It is obvious that 



A4 = o H(N) = VA^(;c). (92) 

However, now the constraints form a Lie algebra since 

{'y(N),M} = 0, 

{M,M} = 0, (93) 



where diffeomorphism constraints nicely form an ideal. The 
master constraint can be regulated via a point-splitting strat- 
egy S as: 



M' 



d/d/ 



XX J X - y) 



H(x) H(y) 



(94) 



Introducing a partition !P of the 3-manifold S into cells C, we 
have an operator i/^„ acting on spin-scalar-network basis Ts^c 
in 'He via a state-dependent triangulation. 



Hc,a ■ Ts,c = Y ^c(v)H{N): ■ r,,. 



(95) 



veV(a) 



where a denotes the underlying graph of the spin-network 
state Ts, and 



H(N):- Z^ci+Z% 



(96) 



i'(A)=i' 



with 



h: 



e.A 



^ 320(v) 
'^*''' 3inyE(v) 
64 



e(./.y.^)e^^^Tr(^-;^(A)^,,(A)[/i;;,^), ^Vf;.]) 



(inyj^Eiv) 

X Tr(/i„(A)[/i;/(A)' K]hsj(A)[h-jf^^^, i']/i,^(A)[/i;](A), y|Vuf^), 

(97) 



and 



H: 



l6N(v) 



0-'(y) 



^■" 3y\mf' 



(98) 



H. 



4,v 



Z 



2'^N{v) 



3^y^{infE\v) 



7r(v) 



v{h)=v(A')=v{X)=v 

X Tr(T,/i,,(A)[/i;;(^,,i]) 

X KiAAKli/^yiyUl)''^]) 
-1 m7__^1/2i 



X h.A^')Ki^'y^^uO'^^ 



1/21 



X /j...(A')[/^;;(A')' (%)']) 



(99) 
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K 



X /^.«(A')[/';;(A')'(^f/j)'''])' 



^6,>. = 2^^"^^^*^^"^^^^^' 



(100) 



(101) 



H: 



7,1- 



3r2/io(/;i)2£(v) 



e(0)=i' v(A)=v 
X [f/io(0(f,,(A))) - I/i„(0(s.,(A)))] 

X Tr(r,/2,,(A)[/i;;(A),(V^,0'^'] 
X /^ma)[/'^;,\a),(V£/,')'^']). 



(102) 



Note that the family of operators //^ ^ are cylindrically con- 
sistent up to difFeomorphism. So the inductive limit operator 
He is densely defined in 'He by the uniform Rovelli- Smolin 
topology. Moreover, the adjoint operators of H^^a' which are 
also cylindrically consistent up to difFeomorphism, read 



iHcJ ■ Ts.c = Yj Xc(v)(H(NrS ■ T,, 



(103) 



veV{a) 



The inductive limit operator, {Hc)\ of {Hq^)^ is adjoint to 
He- Then we could define master constraint operator M. on 
difFeomorphism invariant states as 

m^Diff) ■ T,.c = „ lim ^ ^DifA- y H^ciHc^ ' ^.,.1 

P->S,e,£'->0 2 '^-^ 

ceP 

(104) 

Note that our construction of A1 is qualitatively similar to that 
in 1I27I1 . alt houg h the quantitative actions are different. Similar 
to those in [27] we can prove the following properties of M. 
(i) yV( is diffeomorphism invariant, i.e., 

(f/^AfOay/) ■ r,,, = (M^Diff) ■ Ts,c, 

where U' is induced by the unitary operator in 'Kg corre- 
sponding to a finite difFeomorphism transformation (f. 



(ii) For any given difFeomorphism invariant spin-scalar- 
network state r[5(,], the norm ||yV(r[ic]||£),yy is finite. So VV( 
is densely defined in "Hoiff- 

(iii) yV( is a positive and symmetric operator in "Hoiff and 
hence admits a unique self-adjoint Friedrichs extension. 

In conclusion, there exists a positive and self-adjoint oper- 
ator yVf on "Hoiff corresponding to the master constraint ( |9T1 ). 
It is then possible to obtain the physical Hilbert space of fCR) 
gravity by the direct integral decomposition of 'Hoiff with re- 
spect to M. 



VI. CONCLUDING REMARKS 

How to unify quantum mechanics with gravity theory is one 
of the core problems in modern physics. In recent twenty-five 
years, LQG has made considerable progress in quantizing GR 
non-perturbatively and hence become a fascinating candidate 
theory for quantum gravity. This background-independent 
quantization relies on the key observation that classical GR 
can be cast into the connection-dynamical formalism with the 
structure group of S U(2). Due to this particular formalism, 
LQG was generally considered as a quantization scheme that 
applies only to GR. This was taken by many researchers to 
be a limitation of the quantization scheme. The fact of being 
of general applicability would therefore be significative for 
the general debate about quantum gravity. Especially, fCR) 
gravity theories have become topical in issues related to dark 
energy in cosmology and non-trivial astronomic tests beyond 
GR. Hence, whether such modified gravity theories could be 
quantized non-perturbatively is itself an interesting question. 

The main results of Ref.[10] and the current paper can be 
summarized as follows, (i) The connection dynamics of fCR) 
gravity has been obtained by canonical transformations from 
it's geometric dynamics, (ii) Based on the ,s'M(2)-connection 
dynamical formalism, the rigorous kinematical framework of 
LQG has been successfully extended to metric fCR) gravity 
theories by coupling with a polymer-like scalar field. The im- 
portant physical result that both the area and the volume are 
discrete at quantum kinematical level remains valid for fCR) 
gravity, (iii) While the Hamiltonian constraint is promoted 
to well-defined operator in the kinematical Hilbert space, the 
master constraint can be promoted to well-defined operator 
in the difFeomorphism invariant Hilbert space of loop quan- 
tum fCR) gravity. Thus, the non-perturbative loop quantiza- 
tion procedure is not only valid for GR but also valid for a 
rather general class of 4-dimensional metric theories of grav- 
ity. Therefore, the achievements which have been obtained 
are in two fold. First, classical metric fCR) theories have 
been successfully quantized non-perturbatively. This guaran- 
tees the existence of fCR) theories of gravity at fundamental 
quantum level. Secondly, the valid range of LQG has been 
considerably enlarged to include a rather general class of met- 
ric theories. 

It should be noticed that classically the scalar field (p charac- 
terize different /(■??) theories of gravity hy (p - f'CR). Thus for 
a given fCR) theory, (p becomes a particular function of scalar 
curvature K while the potential ^(<p) is fixed. Hence our quan- 
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turn fCR) gravity may be understood as a class of quantum 
theories representing different choices of the function /(^). 
Of course, there are still many aspects of the connection for- 
malism and loop quantization of /("R) theories which deserve 
discovering. For examples, it is still desirable to find an action 
for the connection dynamics of /(^) gravity. The semiclassi- 
cal analysis of loop quantum fCR) theories is yet to be done. 
To further explore the physical contents of the loop quantum 
/(■??) gravity, we would like to study its applications to cos- 
mology and black holes in future works. Moreover, It is also 
desirable to quantize /(^) theories by covariant spin foam ap- 
proach. 



Since only FJJ^ contains the derivative ofE^ in above equation, 
we consider 



N 



2 

N 






2 h h 



(108) 
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Appendix 

We use {K[j,E^.) and (<^,7r) as canonical variables to derive 
the constraints algebra. By the first canonical transformation. 
The Hamiltonian constraint (fT2b in section 11 can be written as 

20 3yfh <P 



Therefore, we use {k'a{x),E''i^{y)} - 5''g5[6{x,y) to calculate 



{ hvV^(5,0)5,(^), -—{-{k'^eIkIe"^ 

- ^KEfKiE'j) + ^K',Efn)\^^E) -M^N 



+ -^(^KEfn+ln^cl,)+VhD,DPg 
2Sh 3 3 



-L 



and 



(105) 



To calculate the Poisson bracket between two smeared Hamil- 
tonian constraints, we notice that the non-vanishing contribu- 
tions come only form the terms which contain the derivative 
of canonical variables. Those terms are \ cPxN 'sfhDaD"(f), 
which contains both the derivative of Zs^ and the derivative of 
(f>, and J — \(f>N yhR, which only contain the derivative of £*. 
Hence we first use {(l)(x),7r(y)} - 6^(x,y) to calculate 



2£f,2 



■b^iua 



-M{dam{Dc4>)^{-{E''iKlE) 



- ^-ElKiE']) + ^-Eln)) 
EIE', 



+ -M{daN){D,4>)^{-E',)(-{E]Kl'E\ 



^dp 



-E]K'^El) + -E'jn)) -M^N 



(109) 



\Ke'IkIe''j) + p'.Efn)},^^,, -M^N 
- -M{d'N){D,4,)Ke^(-(ElKlE"j 



1 



E'!E'=. 

I I 



r I- r M A -■ 2 . 

{ A^ ^D„D"cl>, -—(-KlEln + -n' 
Jy JY2ih 3 3 

= fiMD,D"N - ND„D"M)(^K<f> + ^KlEf) 



Note also that 



.") 



-M ^ N 



(106) 



4 h 

- \E''jKl'Ei) + ^E''jK))-M^N. 



'h 5-m irrf 



{-3E'^)(-{E]K"^E\ 



(110) 



The combination of above two Poisson brackets equals to 



1 



(107) 



I {NLfM - MD^NX-^nDaCf 
Jn 3 

^{KiE)h„,D>'^-^-KlE)D,^)). 



(Ill) 
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The variation of the terms containing a derivative in C - Taking account of Eqs. (ll06l l-( ll 13b . we obtain 
\(I)N yfhR reads 






'\(I)N) + h"''D,D'{(l)N))6hab 



Dbi^N) - habDcD'i4>N))6li 



ah 






fi yfhiDaDhi(/>N) - habD,D'i,pN))6(^j-^) 



(112) 



Thus we have 






2£f 



= r - ^{MDaDb((f>N) - habMD,D'(cf>N))-j^ 

- \{-2MD,D\(pN)){-E[AE',KiE) 
1 4 

= f - (MDaDb((f>N) - habM{D,D'(pN))h'''\KiE" 



if' 




I - M{DaD''^N)-klE° + M{DcD' 
(—KiEi --k)-M <-^ N 



(f>N) 



= jiNDaD'' 



1 - 



Jpd-^ 



"(0M) - MDaD\cl)N))-KlE) 

+ (ND,D'(<f>M) - MD,D'((/>N))(^n - ^^E". 

= fiNDaD'M - MD,D'N)(^nif> - ^KiEp 

+ (ND,M - MD,N)iD'cf>)djT - ^KiEp 
+ {NDaD'^M - MDaD''N)KlE"j 

2D''d) ~ : „ 

+ {NDaM-MDaN) -K'iE"j. 



(113) 



I {ND.D'M - MD,D'N)(-KiE"i) 

+ (ND"M - MD"N){nD„(p) 
+ iND„D''M - MDaD''N)klE''j 

= UnD"M - MD"N)(D„(KiE'j) - DbiKiE''^ + nD„<f 

- ((DaM)D''N - {D''M)D,N)kIe''j 



L 



(iVg-M-«D°iV)V,-"^"°-"-f'°""''g, (114) 



where we used the following identity 



-((D,M)D^N - {D''M)DaN)klE) 
-{(DaM)DcN - {D,M)DaN)h'"E''.kl 

-2(D^aM)(D,]N)^E"jkl 
-2(DaM)(DcN)^-^kiE"' 



-^'V* 



E^.E' 



e'J\D„M)(D,N)^kl"E"'sk„„ 
n 



[E"DaN, E''DbMf 



■@k 



(115) 



Using above result and shift conjugate pair (k'^,E'') to 
(AJ,, £*), we can easily get the Poisson bracket ( l46b between 
the smeared Hamiltonian constraints. 
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